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Abstract 

We consider a minimally coupled scalar field with a monomial potential 
and a perfect fluid in flat FLRW cosmology. We apply local and global dy¬ 
namical systems techniques to a new three-dimensional dynamical systems 
reformulation of the field equations on a compact state space. This leads to 
a visual global description of the solution space and asymptotic behavior. At 
late times we employ averaging techniques to prove statements about how the 
relationship between the equation of state of the fluid and the monomial ex¬ 
ponent of the scalar field affects asymptotic source dominance and asymptotic 
manifest self-similarity breaking. We also situate the ‘attractor’ solution in 
the three-dimensional state space and show that it corresponds to the one¬ 
dimensional unstable center manifold of a de Sitter fixed point, located on 
an unphysical boundary associated with the dynamics at early times. By 
deriving a center manifold expansion we obtain approximate expressions for 
the attractor solution. We subsequently improve the accuracy and range of 
the approximation by means of Pade approximants and compare with the 
slow-roll approximation. 
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1 Introduction 


The present paper investigates general relativistic flat Friedmann-Lemaitre-Robertson- 
Walker (FLRW) models with a minimally coupled scalar field with a monomial po¬ 
tential, U(0) = ^(A (j)) 2n (A > 0, n — 1, 2, 3,...), and a perfect fluid. The perfect 
fluid is assumed to obey a linear equation of state, p m = (y m — 1 )p m , where p m and 
Pm > 0 are the pressure and the energy density, respectively. The adiabatic index 
is assumed to satisfy 0 < 7 m < 2, where = 1 corresponds to dust and 7 m = 4/3 
to radiation. When 7 m = 2/3 the matter term p m can be reinterpreted as — | 3 R, 
where 3 R is the spatial 3-curvature of the open FLRW model, i.e., 7 m = 2/3 leads to 
equations that are the same as those for a scalar field in open FLRW cosmology. The 
case 7 m = 0 corresponds to a matter content described by a cosmological constant, 
i.e., p m = A, while = 2 describes a stiff perfect fluid; both cases are associated 
with significant bifurcations, and we therefore refrain from discussing them. 

The Einstein and matter field equations for these models are given by 


3 H" — |0 2 + — (A0)^ n + P m ~ P<t> P m ’ (la) 

H = (</> 2 + 7mPm) , (lb) 

0 = 4> + 3H<j)+ A 2n 0 2n "\ (lc) 

Pm 3Lf7mPm- (lb) 


ffere an overdot signifies the derivative with respect to synchronous proper time, t: 
H is the ffubble variable, which is given by H — d/a, where a(t) is the cosmological 
scale factor, and throughout we assume an expanding Universe, i.e. H > 0, where 
H is related to the expansion 0 according to H = 9/3. We use (reduced Planck) 
units such that c = 1 = 87 tG, where c is the speed of light and G is the gravitational 
constant (in the inflationary literature the gravitational constant G is often replaced 
by the Planck mass, G = 1/m^/). 


Heuristically eq. (lc) can be viewed as an equation for an anharmonic oscillator 
with a friction term 3 Htp. This suggests that (0, 0) —> (0, 0) toward the future in 
an oscillatory manner, which is indeed correct. This qualitative picture, however, 
does not show how this comes about in a quantitative way, nor how the fluid affects 
the situation via its influence on H. Running the time backwards allows one to 
heuristically interpret 3 H(j) as an energy input, which suggests that the scalar field 
oscillates with increasing amplitude toward the past, but this picture breaks down 
in the limit H —y 00 . Even though this is beyond the Planck regime, this limit is 
also needed in order to describe the classical behavior at early times after the Planck 
epoch. Furthermore, eq. (Id) yields that 


p m = Po(a/a 0 ) 37, 


( 2 ) 


where po an ao are constants, and hence p —> 0 at late times while p —> 00 at early 
times. 
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Note that the above qualitative considerations say nothing about how e.g. 


r = 


P<t> 

Pm 


( 3 ) 


behaves asymptotically, i.e., if the model is fluid or scalar field dominated, or neither, 
asymptotically. Nor does the above say anything about the role of the so-called 
attractor solution in a global solution space setting. 

This exemplifies that there is a need for a more careful examination, which is illus¬ 
trated by some previous heuristic considerations for a scalar held with a monomial 
potential by e.g. Turner [I] and Mukhanov [2j p. 242, which in turn inspired the 
rigorous work by Rendall j3j; in addition de la Macorra and Piccinclli introduced 
a new heuristic approach to study dynamics at late times for a scalar held with a 
monomial potential and a perfect fluid [ 3 ]; rigorous work in this context was also 
obtained for the special case n — 1 by Giarnbo and Miritzis 0 Q Nevertheless, this 
still leaves room for improvements and extensions, and, as will be shown in this 
paper, it is possible to shed light on interesting previously neglected physical and 
mathematical aspects. 

The main purpose of this paper in, primarily, mathematical cosmology is two-fold: 
Firstly, to obtain a global visual picture of the solutions space, thus, e.g., situat¬ 
ing the so-called attractor solution in a global solution space context. Secondly, to 
prove issues concerning asymptotical behavior at late and early times. This includes 
introducing averaging techniques to determine late time behavior, generalizing and 
simplifying earlier proofs in the literature, and using center manifold theory to rig¬ 
orously derive approximations for the attractor solution at early times, as well as 
clarifying the physically important issue of asymptotic self-similarity. 

The outline of the paper is as follows. In the next section we introduce our new 
three-dimensional dynamical systems reformulation of the field equations on a rel¬ 
atively compact state space. We also present two other complementary dynamical 
systems formulations of the field equations, which allow us to effectively obtain 
approximations for the attractor solution. In Section [3] we apply global and local 
dynamical systems techniques to obtain a complete and illustrative picture of the 
solution space and its properties, including asymptotics. In particular, we introduce 
averaging techniques in our global dynamical systems setting, which allows us to 
prove the following theorem: 

Theorem 1.1. 

(i) If 7 m — > 0, then r = p^/p m —» oo for all solutions with p$p m > 0, which 

implies that the solutions are future asymptotically scalar field dominated. 

(ii) If 7 m — < 0, then r = p<p/p m 0 for all solutions with p^p m > 0, and thus 

the solutions in this case are future asymptotically perfect fluid dominated. 


^ome further examples of references that describe minimally coupled scalar field cosmology in 
dynamical systems settings are Eiia ei, with additional references therein. 
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(iii) If 7 m — = 0, then r = p<j,/p m —» const, when p^p m > 0, and there is thus 

no future scalar held or perfect fluid dominance. 

It should be pointed out that similar conclusions have been reached heuristically 
with quite different arguments in e.g. [4], Furthermore, Giarnbo and Miritzis gave 
a proof for n = 1 for the cases (i) and (ii) in [5] (in the case of general relativ¬ 
ity). However, apart from that our proof rigorously generalizes previous results, our 
method can, in principle, be modified to treat even more general situations. More¬ 
over, we tie our results to the global dynamical systems picture and discuss their 
physical implications, e.g., situating them in the context of future manifest asymp¬ 
totic self-similarity breaking. In Section [4] we focus on the attractor solution, where 
we introduce and compare several approximation schemes, such as center manifold 
and slow-roll based expansions and Pade approximants, in order to describe it quan¬ 
titatively. Finally, Section [5] contains some general remarks, e.g. about the de Sitter 
solution on the unphysical boundary of the state space. 


2 Dynamical systems formulations 


2.1 Global dynamical systems 


Our main global (i.e. compact) dynamical systems formulation is based on the 
dependent variables T,X, Ep which are defined as follows: 


A 0 


(T, X, E t ) , c + R1/n , ( 6nH2 y/2n ’ 

(H, 0, 0) = ( c n f ~ n , V&T~ l X, VQc n T~ n Y^ , 


where 


T = 


T 


1 — T’ 


In addition it is of interest to define 





P<t> 


= E? + X 2n , 


Q r 


Pm 

3 IP' 


(4a) 

(4b) 


(4c) 


( 5 ) 


To introduce a new suitable time variable we note the following: At early times it 
is natural to use a Hubble-normalized time variable r defined by dr/dt = H , due to 
that the expansion 6 = 3 H provides a natural variable scale when 6 —y oo via the 
Raychaudhuri equation, as further discussed in e.g. |9], and references therein (in 
an inflationary context r is often interpreted as the number of e-folds N). At late 
times the square root of the second derivative of the potential, 0 n_1 (for simplicity 
we here incorporate A into 0), provides a natural variable (mass) scale. Due to the 
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Gauss constraint (la), which relates a scale given in 0 to one given in H according 
to 0 ~ H 1 / 71 , this scale can be expressed in terms H according to which leads 

to a dimensionless time variable f defined by df/dt = constant ■ Z/ 1-1 /”, where the 
constant have the same dimension as H 1//n . To incorporate these features in a global 
dynamical systems setting we introduce a new time variable f that interpolates 
between these two regimes, 

f^d-rr 1 , (6) 


where dr/dt —> H, cH 1 1 / n when r —> —oo and r —> +oo, respectively. 

The above leads to the following three-dimensional dynamical system for (T, Ep 


f = ir(l-mi + 9 ), (7a) 

/'7 V' 

—t = -(1 - T)(2 - g)E t - nTX 2n -\ (7b) 

CL I 

dX 1 

— = -(l-T)(l + g)X + TEp (7c) 

dr n 


where the deceleration parameter, q, defined via H = — (1 + q)H 2 , is given by 


Q IT T 7mflm) 1 T 3Sj. T 


( 8 ) 


where 


= Ej + X 


2 n 


^m. 1 


V 2 


X 2n = 1 - % > 0, 


( 9 ) 


where the last equation follows from the Gauss constraint (la), while the inequality 
is due to p m > 0. Above we have also introduced an effective equation of state 
parameter for the scalar field which is defined according to 


70 


1 + P± 
P<j> 


T 2n(^) 2n 


¥ 2 + 27( A< ^) 2n 


( 10 ) 


From the above relations it follows that — 1 < q < 2. In addition it is of interest to 
give the following auxiliary evolution equation for 


^ = 3(1 - T)( 7m - 70)^™, (11) 


2 The variable has been used ubiquitously in the scalar field literature (often denoted by x), 
while X was used in m where it was denoted by y 1 however, as far as we know, the variable 
T and the independent variable f are new. The reason for the name Ej is that mathematically 
this variable plays a role that is reminiscent to that of Hubble-normalized shear, which is usually 
denominated by S in anisotropic cosmology, for a number of situations (the subscript f follows 
the notation in El)- Thus the present nomenclature is designed to pave the way for eventually 
situating the present problem in a broader context than isotropic scalar field cosmology. 
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where 

= 2Y 2 f . ( 12 ) 

The state space S associated with eq. ([7|) is given by a finite (when n > 1 deformed) 
cylinder described by the invariant pure scalar field boundary subset, Q m = 0 (i.e. 
Pm = 0), and thus 0.$ = 1, which we denote by S|n m= o, and 0 < T < 1. From 
now on, we analytically extend S to the state space S by including the unphysical 
invariant submanifold boundaries T = 0 and T — 1. Although these boundaries are 
unphysical, we stress that it is essential to include them since they describe the past 
and future asymptotic states, respectively, of all physical solutions. 

Note that = 0 (i.e. Ef = 0 = A") and hence Q m = 1 is an interior invariant 
subset, S|n o= o, which is just a straight line in the center of the cylinder, describing 
the flat FLRW perfect fluid model without a scalar field (this solution appears as a 
straight line in the present state space due to that it is a self-similar solution, where 
T describes the temporal change in the dimensional variable H). Note also that the 
dynamical system (JT]) is invariant under the discrete symmetry (A", Ef) —>■ —(A, Ef), 
leading to a double representation of the physical solutions when Ul m > 0, which is 
a consequence of that the potential is invariant when (f) —> — 

To describe the dynamics on the scalar field boundary S|n m =o, where = 1, it 
is useful to introduce a complementary global formulation, which is based on the 
following transformation of E^ and X: 

R — ms 2n f) 

E t = F(0)sin0, A" = cos 6*, F(9) = \- -(13) 

V 1 — cos^ a 

This leads to the following regular unconstrained two-dimensional dynamical system: 


^ =-T(l-T) 2 (l-cos 2n £), (14a) 

dr n 

df) Q 

= -TF{6) -(1 - T)F 2 (9) sin 2 9. (14b) 

dr 2 n 

In this case the deceleration parameter q is given by 


q = 2 - 3cos 2 ”9. 


(15) 


The system (14) constitutes a generalization of the system used in [12] . Note that 
for n > 1 the present 9 variable is not the same 9 as that in [3|, which in turn was 
based on the variables used in id. 


3 The system ([7] is differentiable for non-integer n when n > 1, where the differentiability 
depends on n, describing problems with potentials V = 5^(A|</>|) 2r \ where X is to be replaced with 

l*| ^ ®- 
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2.2 Complementary non-bounded dynamical systems 


We here introduce two complementary dynamical systems on unbounded state spaces 
that are useful for describing the dynamics at early times. The first system is based 
on the dependent variables T, A", Ey and the independent variable r, where we recall 
that T and r are defined by 

T = V — = cH~ 1/n % = H, (16) 

1 - T dt v ’ 

where r can be viewed as the number of e-folds N, i.e., r — N. This leads to the 
dynamical system: 

dT 1 ~ 

— = —T(l + g), (17a) 

dr n 

— 1 = -(2 - g)E t - nTX‘ 2n ~ 1 , (17b) 

LL I 

dX 1 

^ = -(l + 9 )V + fE, (17c) 

where q is still given by (|8j) , (J9]) It is also useful to consider auxiliary equations for 
Ll(p and r P<pjPm 


dr 

dr 

dr 



(18a) 

3(7rn - 7 

(18b) 


The second complementary dynamical system concerns the dynamics on the scalar 
held boundary D m = 0. Expressed in terms of T and r the unconstrained system (14) 
takes the form: 


dT 3 ~ 
— = -T 1 
dr n 


-2n 


COS 


0 ), 


df) ^ 

— = —TF(9) - F 2 (0) sin 29. 

dr v ; 2 n V ’ 


(19a) 

(19b) 


Note that the above systems share the same equations as the previous ones with 
bounded state spaces on the invariant boundary subset T = 0 = T. 


3 Global dynamical systems analysis 


It follows from ([7| that 

= 6'nT 3 (l-T) 2 , 

l+q=0 

4 In the special case H m = 0 and n = 2 this system coincides with eq. (16) in [14] : incidentally, 
this model was also the example discussed by Linde in his paper “Chaotic inflation” m- 



dT 

df 


= 0, 


d 2 T 


1+o=0 


dr 2 


= 0, 


1+o=0 


(DT 
df 3 
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on S (note that since we have assumed that 7 m > 0 it follows from eq. & that 


q+ 1 = 0 only when O m = 0 and 7^ = 0). Due to that q + 1 > 0, it follows from (7a) 
and (20) that T is a monotonically increasing function on S (although eq. (20) shows 
that solutions on the scalar held boundary have inflection points when 1 + q = 0) and 
hence T can be viewed as a time variable if one is so inclined. As a consequence all 
orbits (i.e. solution trajectories) in S originate from the invariant subset boundary 
T — 0, which is associated with the asymptotic (classical) initial state, and end at 
the invariant subset boundary T — 1, which corresponds to the asymptotic future, 
and therefore all fixed points are located at T = 0 and T = lj^] It also follows from 
the monotonicity of T that the past (future) attractor resides on T — 0 (T — 1). 

The equations on the subset T = 0 (or, equivalently T = 0) are given by 


<ffi t 

df 

dX 

df 


-(2 -g)S t , 


— (1 + q)X , 
n 


(21a) 

(21b) 


as follows from (JtJ) , or, equivalently 0 (with r replaced with r), where q is given 
by (§ and @. It follows that the state space on T = 0 is divided into four sectors 
defined by the invariant subsets E-j- = 0 and A" = 0. The intersection of these subsets 
with D m = 0 and with each other yield five fixed points on T = 0: 


M±: 

A = 0, 


= ±1, 

(22a) 

dS±: 

A = ±1, 

Ef 

= 0, 

(22b) 

FL 0 : 

A = 0, 


= 0, 

(22c) 


where M± are two equivalent fixed points for which Q m = 0 and q — 2 (and 7^ = 2), 
i.e., they are associated with a massless scalar held state, while the two equivalent 
fixed points dS±, for which Ll m = 0 and q — — 1 (and 7^ = 0) correspond to 
a de Sitter state j^] The fixed point FL 0 gives Ll m = 1 and q = ^(37 m — 2) and 
corresponds to the hat perfect fluid Friedman model. 

As shown below, the two hxed points M± are sources on S; the hxed points dS± 
are sinks on T = 0, but they also have one zero eigenvalue that corresponds to a 
one-dimensional unstable center submanifold on S|n m =o> be., one solution, called an 
attractor solution, originates from each hxed point dS± into S on the pure scalar 
held boundary subset Ll m = 0. The hxed point FL 0 is a saddle that gives rise to 
a 1-parameter set of solutions entering S (the associated unstable tangent space is 

5 A fixed point, sometimes called an equilibrium, critical, or stationary point, is a point xq in 
the state space of a dynamical system x = f(x) for which f(x 0 ) = 0. 

6 Note that the present de Sitter fixed points are distinct from de Sitter states that are associated 
with potentials that admit situations for which dV/d(f> = 0 for some constant finite value of cj> for 
which both V and H have constant, bounded, and positive values. In contrast the present de 
Sitter states correspond to the limits 7 = 0, (0, V, H) —► (±00, 00, 00), and therefore reside on the 
unphysical boundary T = 0. 
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given by = 0), one being the perfect fluid solution given by 
system (22) admits the following conserved quantity when Vt m = 


= 0, Vt m = 1. The 

1 - £2 _ x 2n > 0: 


x (2- lm )n q- 1 = const ) 


which determines the solution trajectories on the T = 0 subset, see Figure [T] 


(23) 



(a) V(4>) = \m 2 (t) 2 , 7m = | 


(b) V(<t>) = j(Ab) 4 , 7m = 1 



Figure 1: The invariant T = 0 boundary subset of S for monomial potentials and a 
perfect fluid with a linear equation of state. The flows are topologically equivalent 
for all monomial potentials with n > 1 and perfect fluids with 0 < 7 m < 2. Depicted 
are two examples. 


The equations on the subset T = 1 are given by 


= nX 2n ~\ 
dr 

(24a) 

dX 

(24b) 

dr * 

This system has a non-hyperbolic fixed point Q 


FLy X = 0, S-j- = 0, 

(25) 


with three zero eigenvalues. Fortunately FLi resides at the intersection of two 
invariant subsets: the invariant subset T = 1 and the invariant perfect fluid subset 
= 0, Q m = 1, and thus, since T is monotone in S, FLi attracts at least this orbit. 
On T = 1 FLi is conveniently analyzed by considering eq. © on T = 1, which 
yields that 

(26) 


_ v) + A = const, 


which is an integral of (24), i.e., the subset T — 1 is foliated with periodic orbits 
surrounding the fixed point FLi, see Figure [2j 


7 A fixed point is hyperbolic if the linearization of the dynamical system at the fixed point is a 
matrix that possesses eigenvalues with non-vanishing real parts; if the linearization leads to one or 
more eigenvalues with vanishing real parts it is said to be non-hyperbolic. 
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1% ret 




Figure 2: The invariant T = 1 boundary subset of S for two examples: The poten¬ 
tials V(cp) = and V((p) = |(A0) 4 and a perfect fluid (with arbitrary equation 

of state 0 < 7 m < 2). 

To make further progress as regards the global properties of the solution space we 
need to consider the asymptotic dynamics at early and late times. 


3.1 Asymptotic dynamics at early times 


The linearization of the system ([7]) at the fixed points FLo and My is conveniently 
described as follows: 


1 dT 
T df 
1 dSj 
Sj df 
1 dX 
X ~df 


FL 0 


FL 0 


FL 0 


3 

1 

dT 

7 7m) 

2n 

T 

df 

— 2 _ 7m) 

1 d( 1 ± 

1 ± Ej. df 

E t ) 

3 

1 

dX 

7 7 mi 

2n 

X 

df 


M± 


Mn 


M± 


3 

n 


3(2 - 7m), 


3 

1 

n 


(27a) 

(27b) 

(27c) 


where the right hand sides constitute the eigenvalues of the fixed points. Hence FL 0 
is a hyperbolic saddle, with an unstable manifold tangential to = 0, i.e., there 
is a 1-parameter set of solutions that originate from FL 0 entering the state space S 
tangentially to E| = 0. The fixed points M± are hyperbolic sources; it follows from 
the invariant submanifold structures that there exists a 2-parameter set of solutions 
entering the interior of the cylinder with O m > 0 from each fixed point M±, while a 
1-parameter set of solutions originate from each fixed point M± into the boundary 
subset S| Qm=0 Q 

s The above results follow from the Hartman-Grobman theorem, which states that in a neigh¬ 
borhood of a hyperbolic fixed point the full nonlinear dynamical system and the linearized system 
are topologically equivalent, see e.g. mm- 
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On the subset T — 0 linearization of the system (22) gives 

1 dEj 


dr 

d( 1 =f X) 


— — 3 , 


dS± 


ItX 


dr 




dS± 


(28a) 

(28b) 


and hence dS± are hyperbolic sinks on T — 0, as illustrated in Figure [lj In the 
full state space, however, each equivalent fixed point dS± have an additional zero 
eigenvalue associated with a one-dimensional so-called center manifold. Fortunately, 
the center direction lies on the S|f2 m= o subset, and hence we can investigate the 
center manifold by means of the unconstrained system (14), which will be done 
in Section [4j There we show that the center manifold of each (equivalent) fixed 
point dS± corresponds to a single solution that enters the state space S|n m=0 (we 
will even obtain approximate expressions for this solution), and this solution, which 
hence resides on S|o m= o, is what is often referred to as the ‘attractor’ solution. In 
the full state space S the fixed points dS± are thus center-saddles. 


From these considerations, in combination with the monotonicity of T, it follows 
that all solutions are past asymptotically self-similar in the sense that all physical 
geometrical scale-invariant observables, such as the deceleration parameter q, are 
asymptotically constant. However, there is a twist to this. The geometry of a flat 
FLRW pure perfect fluid cosmology with a linear equation of state and the geom¬ 
etry of a pure massless scalar field are geometries that admit a proper spacetime 
transitive homothety group, and such spacetimes are invariant under scalings of the 
spacetime coordinates. This is the underlying reason why they can be represented 
by fixed points, but not all fixed points are associated with geometries admitting 
proper spacetime transitive homothety groups, as exemplified by the de Sitter space- 
time. For these spacetimes homothetic scale-invariance is broken by the dimensional 
cosmological constant, but the 1-parameter set of de Sitter spacetimes (parameter¬ 
ized by A) admits a scaling (self-similar) property that scales A (i.e., a scaling that 
maps one de Sitter spacetime to another with a different A), and it is due to this 
scaling property de Sitter spacetimes can appear as fixed points. 


In the present case the fixed points are not in the interior physical state space, but 
on the unphysical boundary, but they are nevertheless characterized by e.g. the 
same value of q as the associated physical spacetime. The fact that they in the 
present context appear on the unphysical boundary has consequences, which we will 
come back to in a discussion about the de Sitter fixed points in Section [5] Finally 
we point out that, apart from a set of measure zero, all solutions originate from a 
massless scalar field state and hence the present models are past generically massless 
scalar field dominated. 
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3.2 Asymptotic dynamics at late times 


The equations on the T = 1 subset, i.e. equation (24), are equivalent to that of (lc) 


when setting H = 0, i.e. this problem is exactly that of an anharmonic oscillator 
(when n > 1; for n — 1 the problem is that of a harmonic oscillator). This can be 


seen from eq. (24), which yields 


d 2 X 

df 2 


nX 2 "' 1 = 0. 


(29) 


We will now apply the approximate ideas in Mukhanov [2] to the present exact 
problem of an anharmonic oscillator. We therefore first multiply the above equation 
with X and rewrite it as 


- (x—) 

df\ df J 


£V + n *- = o. 

cir J 


(30) 


Each periodic orbit is characterized by a constant value of and has an associated 
time period P = P(f^). The time average of a function / over a period for a 
periodic orbit characterized by is given by 




0T0 + -P(W) 


fdr/Pff1$). 


(31) 


'TO 


Taking the time average of eq. (30) for a periodic orbit gives 




(32) 


where we for notational convenience from now on drop the subscript Again, 
note that in contrast to the result in [2], this is an exact relation on the subset 
T — 1. Using this result on T = 1 for a periodic orbit in the definition of 7^ yields 


(7</>) = 


\ Ll<f> 


2(Zf) 2<E?> 


2<£?> 


O 


* (EJ + X 2 ”) (E 2 ) + (X 2 ”)’ 


(33) 


which together with (32) leads to 


(70) = 


2 n 

n + 1 


(34) 


on the subset T — 1, i.e., (7^) n = (7is independent of 0$. It therefore follows 
that on average, in the above sense, e.g. n — 1 and n — 2 on T =1 correspond 


to dust and radiation, respectively. Note that the result (34) coincides with the 


approximate heuristic results using proper time given in [T] and [23; see also [3J for 
a quite different precise definition of (7^), which still yields 
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Before continuing it is instructive to consider a model that consists of two perfect 
fluids with constant equation of state parameters 71 and 72 . Then r m = pi/p 2 oc 
a 3 ( 72 — 71 ) (this expression follows from that p oc a -37 , but it can also be obtained 
from the equation dr m /dr = 3(72 — 7 i)r m ). Since it is not difficult to show that 
a —* 00 toward the future it follows that r m —>• 00 if 72 > 71 ; r m —* 0 if 72 < 71 ; 
r m —> const, if 72 = 71 , i.e., the fluid with the softest equation of state dominates 
toward the future. 


Assuming that asymptotically 7 ^ can be replaced with the asymptotic averaged 
result (y^) = in eq. (18b) results in 


dr 

dr 


3 (7m - (7 4>)) r 



2 n \ 
n + l) r ’ 


(35) 


which suggests Theorem 1.1 given in the Introduction, but this is of course no proof. 


Next we introduce averaging techniques that are subsequently used to provide the 
proof of this theorem. 


Averaging 


Standard averaging techniques and theorems can be found in Chapter 4 in (33] (the 
periodic case) and in [32] (the general case). In standard averaging theory, a per¬ 
turbation parameter e plays the key role: roughly speaking, a differential equation 
of the form x = £f(x,t,e) for e > 0 is approximated by the averaged equation at 
e = 0, i.e., y = s{f(y,., 0)), where the average (.) is defined in eq. (31). Further¬ 
more, the error \x — y\ has to be controlled. In the problem at hand, we consider 
the differential equations in the variables fand T, where the role of the parameter 
£ is played by 1 — T. Therefore, after setting 


£ = 1 - T, 


(36) 


we have to prove an averaging theorem for the case where £ is not a parameter, but 
a variable that slowly goes to zero. The evolution equation of 0.$ is given in eq. (18), 
which in terms of £ takes the form 


dO,p 

df 


3£(7m 7</>)^</>(l ^</>)> 


(37) 


where 7 ^ = 2'E^/Ll ( j > . This formulation is problematic due to that 7 ^ is not well- 
defined when Lltj) is zero. We therefore use the following formulation: 


7W 3«r( 7m U> - 2EJ)(1 - SV), 
£ = -!**(!- e)( l + 9) , 


(38a) 

(38b) 
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where (X, solves the system ([7]), + X 2n , and 

1 + 1 — 2 ( 2 ^t + 7m(l — ^)) • (39) 

The general idea of averaging is to express as 

0^ = y + ew(y,e,f), (40) 


and prove that the evolution of the variable y is approximated at first order by the 
solution y of the averaged equation. For that, consider the average as defined in 
eq. (31) of the right hand side of eq. (38a). More precisely, considering an equation 
of the form y' = ef(y,e,f) + &P 2 ) with f-periodic /(y, 0 ,r) of period P = P(y), 
the averaged equation is given by y' = (/)(?/), where (f)(y) p J 0 f(y,0,f)df. 
According to (34), we have 2(Y 2 ) = ( 7 ^)^, where ( 7 ^) = is a constant that 
does not depend on f 1^. Hence the averaged equation reads 


Jr = 3^(7 m - (70»y(i - y ), 


(41) 


while w in eq. (40) will be chosen appropriately in the proof. 


Proof of Theorem 11.II 


To prove Theorem 1.1 we first re-express the theorem in terms of 


(i) If > 0, initial conditions with positive < 1 and e converge 

for t —> 00 to the outer periodic orbit with = 1 tangentially to the slice 

= o}. 

(ii) K 7 m - < 0 , initial conditions with positive 1 — 0,$ and e converge for 

t —* 00 to the center with Vt^ = 0 tangentially to the slice {e = 0}. 

(iii) If 7 m — = 0, each periodic orbit on the slice {e = 0} attracts a 1-parameter 

set of trajectories with positive initial e. 


Proof. Let us first derive a differential equation for y by taking the time derivative 
of eq. (|4C)]): 


df 


dy de 


f dw dy dw dw de \ 


dr dr W £ \dy dr dr de dr) 
dw\ dy dw de dw de 

1 + £ ay)sF + £ aF + *“' + £ &*' 


On the other hand, 
(AH 


dr 


3 e (7 m - <7 0 ) + (70) - 7 <a) (y + ew){l -y- ew) 

3 ^( 7 m ~ < 7 < 7 >) 2/( 1 ~y) + 3 « 7 <i>)y ~ 2£ 2 )(1 - y) 

+ 3e 2 (7m - - 2 y) - 3£ 3 (7m - 7 P)w 2 . 


(42) 


(43) 
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Let us now set 


»r = 3«7*>y-2E?)(l 


df 
^2 


y)- 


(44) 


Note that for large times 2E 2 is well approximated by periodic functions with an 
average (7^)?/. The right hand side of (44) is for large times almost periodic and 
has an average that is zero so that the variable w is bounded. 


As a consequence, the differential equation for the variable y takes the form 


dy 


dr 


dw\ 


-i 


TV — 1 + £tT“ 


dy J 


|3e(7m - <7*»2/(l ~y) 

+ £ 2 W ^3(q m - 7</,)(l - 2 y) + ^(1 - e)(l + g)^ 


(45) 


dw de 
de df 


3e 3 (7 m 



Using the fact that = 0(e 3 ) and that ^1 + e^j = 1 — £§7 + 0(e 2 ) results 

in the following: 


— = 3£(7 m - (7^))2/(l - y) 

f 1 dvo 'l 

+ £ 2 |3w(7 m - 7 0 )(1 - 2 y) + -(1 - £)(1 + q) - 3—( 7 m - (7^))y(l -y)\ 

+ 0{e 3 ). 

(46) 


Next we have to prove that the solution y of this equation and the solution y of 
the averaged equation (41) have the same asymptotics when f —> +00. Since the 
averaged equation (41) is expected to govern the dynamics, we first study the late 
time behavior of the system 


— = 3£(y m - (70))y(l - y ), 
* = _U ( i_ £)( i +9) . 


(47a) 

(47b) 


After Euler multiplication by l/£ (or equivalently, a singular change of time variable 
ed/dr = d/df ), this system reads 


^ = 3(7 m - (7*))j/(l -y), 
de l 

£ = ~n £{l ~ e){l + q) - 


(48a) 

(48b) 
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In cases (i) and (ii), for which 7 m — ( 7 ^) 7 ^ 0, the two fixed points of this system are 
located at (e = 0 ,y — 0) and (e = 0 ,y = 1). The line £ = 0 is a heteroclinic orbit 
between these two fixed points, whose direction depends on the sign of 7 m — ( 7 ^). 


Undoing the Euler multiplication does not affect the trajectories with positive e. 
On the other hand, in the original averaged system (47), the line {e = 0} is a line of 
fixed points. Solutions with positive initial £ will approach the line of fixed points 
at £ = 0 , but slowly slide along this line as t —>• 00 in the direction prescribed by 
the sign of — ( 7 ^), and go to the left or right fixed point accordingly. 


Next we prove that the variables y and ffollow this evolution. First note that the 
sequences {£ n }neN> {T n }neN> defined as follows, 


To = 0 , 

£0 > 0 , 


r n+l T n Y 1 /£ni 

&n +1 4 Ai+1 ) 1 


(49) 


have limits 


li m n —>00 T n H“CXD, 

linin—Kx) 0, 


(50) 


since e(f) goes to zero when f goes to infinity. For a sufficiently small £ > 0, 
eq. (46) guaranties that y is monotone (in- or decreasing, depending on the sign of 
the quantity y m — ( 7 ^)) and bounded. Hence y(f) must have a limit when f —* 00 . 
Next we estimate ((f) := y(f) — y(f), where y and y are trajectories with identical 
‘initial’ conditions at time r n : 


ICWI = | / 3 £( 7 m - < 70 >)(y -y)(l-(y + y)) ds 


-V- 

!l<i 


+ j (e 2 h(y,w,£,s)+0(e 3 )) ds | 
l-l <M 

FT pf 

<3 Ce n \((s)\ds + e 2 n Mds + 0(e 3 n )) 


(51) 


< 3 Ce n / \((s)\ds + £ 2 n M(r - r n ) + 0(4), 


for r > r n+ 1 , where we used the fact that |( 7 m — < 7 ^ >)| < C with C > 0 constant. 
Applying Gronwall’s Lemma (see Lemma 4.1.2 in [23], or p. 37 in Til), results in 


KW | < y£ n (exp(3C'£ rl (T - T n )) - 1) + 0(£ 2 ). (52) 

Hence for f—f n G [0, l/£ n ], i.e. t G [f n , f n+ 1 ], the inequality |^(r)| < K£ n holds for a 
positive constant K. Letting n go to infinity implies that ( tends to zero when f goes 
to infinity. Therefore y and y have the same limit as r — > oo, i.e. 0 or 1 depending 
on the sign of the quantity — (y^). Finally, recall that = y + £w ; from the 






3 GLOBAL DYNAMICAL SYSTEMS ANALYSIS 


17 


triangle inequality, and £ —» 0 when t —> oo, it follows that also converges to 0 
or 1 according to the sign of — (7^). This completes the proof of the theorem for 
the non-critical cases (i) and (ii) for which 7 m — (7^) 7^ 0. 

Let us now consider the critical case (iii), for which 7 m = (7^) = 2 n/(n + 1). In this 
case the right hand side of the averaged equation (41) vanishes. Furthermore, the 
evolution equation (46) for y becomes 

dy 


dr 


= £ 


— 7</>)(l — 2 y) + —(1 — e)(l + q) \ + 0(e ) 


(53) 


Let us first consider the average of the right hand side. We therefore define 
g(y,w,e,f) = 3w('y m - 7^(1 - 2 y) + -(1 - e)(l + q), 


(54) 


and compute its average at e = 0, where £ 2 is a periodic function and (2£ 2 ) = (7 


{g){y,w) = p 9 (y,w, 0 ,r)dr 


= 3w (7 m - (7^))(1 -2y) + -(l + q) 


—v-- 

= 0 


3 1 


= 2 npL (2S t +7m(l ~y))df 


(55) 


= 2 ^ I ~ ^))y + 


= 0 


2n 7m 


Note that in the critical case, ^7 m = ^(7</,) = 2 nj^A) = 7+r We th ere f° re obtain 
the following averaged system, expected to give the leading order approximation: 


dz 2 3 

df n + 1 ’ 
de 3 

df 


71(71 + 1 ) 


£ 2 (l-£). 


(56a) 

(56b) 


Again we study the dynamics on a Euler multiplied version of this system, i.e. in 
the time variable ed/dr = d/df: 


dz 3 
dr n+ 1 ’ 

de 3 


dr n(n + 1) 


e(l-£). 


(57a) 

(57b) 
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The linearization at the line of fixed points {e = 0} is 


0 


O+i) 

3 

n(n+ 1) 


This linearization admits one eigenvalue that is zero with a corresponding eigenvec¬ 
tor that is parallel to the line of fixed points, and a stable (i.e. negative) eigenvalue 
A = 


' n(n+i) i w ith a corresponding eigenvector (z = —n,£ = 1) pointing toward the 
inside of the cylinder (i.e. the line of fixed points of the system (57) is transversally 


hyperbolic). Hence each fixed point (zo,0) has a one-dimensional stable manifold, 


i.e., there exists a trajectory z(r) that solves the averaged system (57) with a positive 
initial e and converges to (zq, 0) for each z 0 as r —>• oo. This asymptotic behavior is 
not affected by transforming the equations back to the time variable f. 

A straightforward estimation of the term 0(e 3 ) provides bootstrapping sequences 
{r n } neN and {£ n } ne N, defined similarly as in the non-critical case (49). In other 


words, we obtain a pseudo-trajectory (O^(r)} of the original system (38) with 


^(r n ) = z{r n ), Vr G [r„,r n+ i], |^(r) - z(r)\ < Ke n , (58) 

for a constant K. By regularity of the flow and compactness of the cylinder, there is 


an initial data whose trajectory H^(r) under the flow of the original equation (38) 
shadows the above pseudo-trajectory in the following sense: 


Vn G N, Vr G [r n ,T n+1 ], |0^(r) - fi</,(r)| < Ke n . (59) 

By the triangle inequality, |f^(r) — A(t)| —)■ 0 as t —* oo; therefore, for each 
Zo G [0,1], there exists a trajectory that is limiting to the periodic trajectory at 
e = 0, characterized by = zq. Translating these results into the state space of ([7]) 
concludes the proof of case (iii). □ 


Expressing these results for orbits with 0 < Lt m < 1 and 0 < T < 1 in terms of 
formal global future attractors A+ of the global dynamical system ([7]) leads tofl 

(i) A+ = S|r=i,n m =o if 7 m > 

(ii) A+ = S| r=1) n,= 0 = FL, if 7m < 

(iii) A + — S\ T =i if 7m = ^y- 

The above results are illustrated by the numerical examples that are depicted in 
Figure [3} 

9 Loosely speaking, in dynamical systems theory attractor behavior describes situations where a 
collection of state space points evolve into a certain ‘attractor’ region from which they never leave. 
For a formal definition of a dynamical systems attractor, see e.g. na. and references therein. 
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Figure 3: Solutions for the scalar field potential V(<f>) = |(A0) 4 for various matter 
equation of states. The first picture describes a solution and the attractor solutions 
on the scalar field boundary. The other figures depict a typical solution that illus¬ 
trates the behavior at late times for cases (i) (i.e., (ii) (i.e., 

and (iii) (i.e., = ^). 

Physical Interpretation 

The present results have physical consequences. The solutions approach (H, </>, (j>, p m ) = 
(0, 0, 0, 0) toward the future, which is the Minkowski solution of the system ([Tj). The 
Minkowski solution admits an 11-dimensional homothety group and is therefore an 
example of a self-similar spacetime. However, the asymptotic future state for solu¬ 
tions in case (i), with and > 0, is the periodic orbit at = 1 — T, 
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and hence, since the future asymptotic behavior is described by a limit cycle, it 
follows that the deceleration parameter q oscillates toward the future (this is also 
obviously true for the orbits on the subset = 1). The solutions thus approach 
the Minkowski spacetime in a foliation that is induced by the present models in a 
manner where the self-similar nature of the Minkowski spacetime is not manifest 
and thus we say that these models exhibit future asymptotic manifest self-similarity 
breaking 

Physically, case (ii), for which < ^kt-, has the simplest asymptotic regime. The 
future state for this case when Ll m > 0 is the fixed point FLi for which q = |(37 m —2). 
This state describes the flat perfect fluid model with equation of state parameter y m . 
Since this model is self-similar, the solutions asymptotically approach the Minkowski 
spacetime in a manifest self-similar manner and we therefore say that the present 
class of models are future asymptotically (manifestly) self-similar. 

Case (iii), where implies that (with the exception of the perfect fluid 

solution with = 0) all solutions oscillate asymptotically toward the future and 
thus this case is also characterized by future asymptotic manifest self-similarity 
breaking. Since asymptotic manifest self-similarity, and breaking thereof, plays an 
important role in physics as a whole, it is of considerable interest to note the above 
features. 

Furthermore, recall that the case = | is equivalent to a pure scalar field case in 
an open FLRW model. As a consequence all monomial potentials (with n > 1) for 
these models lead to that the future end state is given by the fixed point FLi, which 
in this context is associated with a Milne state. 

We are now finally in a position to describe the global solution space and its features. 
There is a 1-parameter set of solutions that originate from the fixed point FL 0 , cor¬ 
responding to an initial perfect fluid dominated state, into the state space S; two 
equivalent 2-parameter sets of solutions that originate from each of the equivalent 
fixed points M±, corresponding to a self-similar initial massless scalar field domi¬ 
nated state, where two equivalent 1-parameter subsets, belonging to each of these 
sets of solution, reside on the scalar field boundary = 1; finally, there are two 
equivalent attractor solutions residing on Qy = 1 that originate from each of the 
fixed points dS-t, thus exhibiting an initial de Sitter state in the sense previously 
discussed. The future state of all orbits in S resides on the subset T — 1, where the 
precise asymptotic location on T = 1 depends on the relation between the mono¬ 
mial potential exponent 2 n and the equation of state parameter according to 

10 Asymptotic (continuous) manifest self-similarity here implies that physical geometrical scale- 
invariant observables such as the deceleration parameter q take asymptotic constant values, but 
since the future attractor is a limit cycle it follows that q is asymptotically oscillating, i.e., asymp¬ 
totic manifest self-similarity is broken. For another example of future asymptotic self-similarity 
breaking in cosmology, see [191 . Note that asymptotic manifest self-similarity is a complicated issue 
in the present asymptotic Minkowski case due that the only metrics with no conformal scalars (and 
hence the only admitting a local conformal group not conformally isometric) are either conformal 
to the plane wave metric with parallel rays or conformally Minkowski, see UBEI]. 
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Theorem 11,11 

In all three cases (i), (ii), and (iii) there exists an open set of solutions that is 
close to the attractor solutions at some intermediate stage of their evolution, but 
there also exists an open set of solutions that is not close to the attractor solutions. 
Furthermore, although there exists an open set of solutions that is close to the 
attractor solutions at late times in case (i) this is not true for case (ii) where each 
(equivalent) attractor solution only acts as a kind of ‘saddle.’ In addition, note that 
the variables T,Ef, X can be locally expressed in terms of the physical quantities 
H, q, Q m . By imposing a Euclidian measure on the space described by II. q. Q rn one 
finds that ‘most’ of the evolution of ‘most’ solutions is not governed by the ‘attractor’ 
solution; the addition of a degree of freedom, in this case associated with a perfect 
fluid, has aggravated the situation for arguments that attempt to establish that 
attractor solutions in some sense are attractors]^] Nevertheless, attractor solutions 
are likely to continue to generate considerable attention, and thus there is a need to 
describe them accurately, which is the topic of the next section. 


4 Attractor solution approximants 

In this section we will introduce and compare several analytical approximation 
schemes that describe the attractor solution. Since the two equivalent attractor 
solutions are just the center manifolds of dS-t, we begin with a direct approach of 
approximately obtaining these center manifolds by means of center manifold theory. 
For presentations of center manifold analysis, see e.g., ESI HE!, and for applications 
in cosmology, e.g., HU [281112]. 


4.1 Center manifold approximants 


The equivalent center manifolds of dS± reside on the scalar field subset Q m = 0 
and therefore it suffices to study the equations on this subset. In [12] it was shown 
that for quadratic potentials the center manifold expansion using the bounded sys¬ 
tem (14) for n = 1 resulted in a larger range than the expansion based on the 
unbounded system (19). However, it was also shown that by using these expansions 
to produce so-called Pade approximants resulted in much better approximations, 
both as regards accuracy and range. Moreover, the Pade approximants for the sys¬ 
tem (19) gave the same expressions as the Pade approximants obtained from (14) 
(or, more accurately, they produced the converging subset of Pade approximants 
associated with (14)), but in a simpler form. Thus the most useful results from 
center manifold analysis are obtained by performing such an analysis for the un- 


11 For further discussions about the meaning of ‘attractor solutions’ and measures, see [T2] and 
the recent papers by Remmen and Carroll [221 S3] and by Corichi and Sloan (24] , and references 
therein. 
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bounded system (19), followed by an introduction of Pade approximants, instead of 
using (14). 


Due to the discrete symmetry, we can, without loss of generality, restrict the center 
manifold analysis to the fixed point dS + at 9 = 0 = T. A linearization of the 


system (19) at this fixed point yields 


E S = {(T,9)\T = 0}, (60a) 

E c = {(T,0)\y/nf + 39 = 0}, (60b) 


where E s and E c denote the tangential stable and center subspaces, respectively. 
Since the tangential center subspace is given by a JnT + 39 = 0, we introduce 

v = T+^=9 (61) 

yjn 


as a new variable that replaces 9 in order to study the center manifold W c , with 
the tangent space E c at (T,v) = (0,0). As follows from (19), this leads to the 
transformed system 


^ =-f (1-cos 2n 0) , (62a) 

dr n v ' 

= cos2 " e) ~ fm - l FH0) sin 20 ) • (62b) 


where 9 = ^(v — T). The linearization of eq. (62b) yields = — 3v, while eq. (62a) 


only has higher order terms. The center manifold W c can be obtained as the graph 
v = h(T) near (T,v) = (0, 0) (i.e., use T as an independent variable), where /r(0) = 0 
(fixed point condition) and y|(0) = 0 (tangency condition). Inserting this relation¬ 


ship into eq. (62) and using T as the independent variable leads to 


~^=T (l — cos 
Jn 


2 n 


«(?-■) 


^ - 1 ) + TF{9) + — F 2 (9) sin 2 9 = 0, 

2 n 


(63) 


where 8 = ^f(h(T) - T). 

Solving the above nonlinear ordinary differential equation amounts to finding the 
attractor solution, which for the present class of problems does not seem likely to be 
possible. Instead the equation can be solved approximately by representing h(T) as 
a formal truncated power series and by Taylor expanding the expressions involving 
9 (and subsequently replace 9 with its power series expression in T), which makes it 
possible to algebraically solve for the coefficients in the formal series. Before doing 
this, it is useful to note that all coefficients in the above equation are odd in terms 
of T and 9 , and as a consequence the power series for h consists of odd powers of T, 
since it is only then ^ results in even powers (odd powers for and hence even 
powers for h, are clearly zero due to the above properties). Furthermore, it follows 
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from the tangency condition that the series for h in odd powers of T begins with a 
cubic term. Hence 

n n 

h(T) = a iT 2i+1 + 0(T 2n+3 ) = TJ2 + 0(T 2n+3 ) as f 0, (64) 

2=1 2=1 


where the series for h(T) is truncated at some chosen order and where we have in¬ 
troduced x = T 2 . Inserting this into a Taylor expansion of eq. (63) and algebraically 
solving for the coefficients leads to 


6 


Tf(x), 


Ti V Ti 

fix) = 1 +- (3 n- 7)x +-(65n 2 - 570n + 1081)x 2 

JK J 108 L V ; 360 V . 


where we have chosen to truncate the series for 6 at 5th order in T. 


(65) 


To improve the range and accuracy of the above approximation for the attractor 
solution we construct the [1 /lU(x) Pade approximant, which leads to the following 
approximate expression for 0\\ 


\ T) ~ 

6 «-V-Ul/l]^), 


[i/i]/W 


1_ n _ (IrR _ 43-n 27531 

36(3n-7) \2 U ^ 30 ) X 

1_ n _ ( Hr;2 _ 37 ?7 1081) T ' 

1 36(3n—7) l 2 Ollier 10 ) X 


( 66 ) 


When this is subsequently expressed in T, i.e., T = T/(l — T) and x = T 2 /(1 — T) 2 , 
this yields a curve 9(T) which approximates the attractor solution in the state space 
S on the boundary Q m = 0. The case n — 1 was dealt with in Figures 5 and 
14 in [CP], where these approximations curves for the attractor solution, as well 
as those obtained by considering higher order expansions, were compared with the 
numerically computed attractor solution. To obtain explicit curves in our state 
space picture we need to specify n. To avoid details, we in this paper only compare 
the approximate solution curves with the numerically computed attractor solution 
for some representative values of n for the above [1/1] f(x) Pade approximant in 
subsection 44 below; although it should be pointed out that higher order Pade 
approximants give better, although more complicated, results. 


Ref. [12] also illustrates that nonlinear transformations might have two important 
features. First, they affect analytical or numerical computations and therefore a 
suitable choice of variables can make a problem more tractable. Second, they affect 
approximations and a suitable choice of variables can lead to better approxima¬ 
tions for solutions like the attractor solution. We will therefore next consider an 
approximation scheme based on the variables Fq and X instead of 6. 

12 A Pade approximant of order (m,n) of a function f(x), denoted by [m/n]f(x), is associated 
with a truncated Taylor series / rj Co + C\X + C 2 X 2 + • • • + c m+n x m+n and given by the polynomials 
P m {x) = Po + PiX + P 2 X 2 + • • • + PmX m and Q n (x) = q 0 + q\X + q 2 x 2 + ■ ■ ■ + q n x n according 
to [m/n]f(x) = , such that Q n (x)(co + c\x + c 2 x 2 + • • • + c m+n x m+n ) = P m (x), where 

coefficients with the same powers of x are equated up through m + n. For details and examples, 
see mmmm and references therein. 
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4.2 Series expansion approximants 


As in the previous case we use a system adapted to the dynamics at early times 
and we therefore use T, but 6 is replaced with Ef and X, i.e., we consider the 
system ©• We are interested in hireling new approximations for the attractor 
solution, and inspired by the result in eq. (65) we assume that Ej and X can be 
written as formal truncated series in T. These series are subsequently inserted 


in (17), but since we are interested in the attractor solution, which resides on the 
fl m = 0 boundary, we also require that the constraint X 2 + X 2n = 1 is satisfied. The 
analysis is simplified by noting that it follows from the system 0 that Xq must 
have only odd terms while A" have only even terms in their series expansions in T. 
Moreover, since we choose, without loss of generality, to consider the solution that 
originates from the dS + fixed point it follows that to lowest order Xj = 0, X — 1. 
Algebraically solving for the coefficients, and writing T 2 = x as before, leads to 


X 


T) ~ 71 71* 

- Tf(x ), f(x) = 1- x-\ -(17 — 6n)x 2 , 

3 JK h Jy ’ 18 648 V ; ’ 

n n 2 . . 2 

1 - Is 1 + 648 (5 - 2n)x ’ 


(67a) 

(67b) 


where we have chosen to truncate the series for f(x) and X(x) at 2nd order in x. 

To improve the range and accuracy we calculate the [1/1]/(x) and [l/l]x(^) Pade 
approximants, which yield the following approximate expressions 

1 + £(5-2 n)x 


n~ 

—T 

3 


1 + £(17 — 6 ri)x 


X 


1 + £ (3 - 2 n) x 
1 + £ (5 — 2 n) x 


(68a) 

(68b) 


Note that for small values of T in the neighborhood of dS + we obtain X = 1— £T 2 = 


cos# ~ 1 — \6 2 , which gives 6 = — ^T, which is just the tangency condition for the 


center manifold of dS + . Replacing T and x with T in (68) yield two curves in the 


state space S on the scalar held subset Ll m = 0, one for the approximation of A" and 
one for the approximation of Ep both approximating the attractor solution. 

To obtain accurate numerical results for the attractor solution it is preferable to use 
the unconstrained system (14) rather than ([7]) subjected to the constraint X 2 +A" 2n = 


1, especially when since the constraint surface then becomes unstable for 

sufficiently large T. The above results are translated into the corresponding 6{T) 
curves via 6 = — arccos X (the minus sign is due to that 6 < 0 and taking the default 
branch) and implicitly via F{6) sin# = —| (yzy) /, where X(x(T)) and f(x(T)) are 
obtained by replacing x with T in (67) and (68), which lead to approximations that 


subsequently can be compared with the numerically computed attractor solution. 
Below, in subsection 4.4| we for brevity only discuss the Pade approximant given in 
eq. (68a) for some representative values of n. 
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As a final remark, note that the lowest order expansion just gives X — 1, which is 
equivalent to = 0, since X 2n = 1 — T, 2 on the scalar field boundary. The condition 
Sj = 0 and X = 1 describes a straight vertical line in S, corresponding to q = — 1, 
originating from dS + , which apart from dS + is not a good approximation for the 
attractor solution. This is in contrast with the previous center manifold expansion, 
which is given by a unique curve d(T), which to all orders in the limit of small T is 
tangent to the center manifold. 

Next we turn to approximations based on the slow-roll approximation and its ex¬ 
tensions. 


4.3 Slow-roll based approximants 


We here extend the work in [S5| and [12] to general n and to higher order by 
approximating i/ 2 (0) with so-called slow-roll Hubble expansions, and illustrate this 
class of approximations in our global state space picture. To facilitate comparison 
with the inflationary literature, we initially keep the coupling constant k = 8nG = 
87rmp 2 in the slow-roll Hubble expansion formulas below. The approach developed 
in |25| is based upon a hierarchy of slow-roll parameters, the first being 


£h 



W \ 

0 2 + V{4)) 



(69) 


which are assumed to be small. These parameters allow one to produce a truncated 
expansion for the Hubble variable in terms of 0, which for the present monomial 
scalar field potential is given by (following the prescription given in [25], which we 
refer to for details): 

3 H 2 » £(A <t>) 2n f{y) = k,V ((/>) f (y), 

2 n 

(70) 


where we have found it convenient to define a quantity y according to 


V = 


2 n 2 

3 K(j) 2 


(71) 


Note that the 2nd (3rd) order term is zero in (70) when n — 2 {n — 3). The 
above series expansion contains no parameters and thus it can only describe a single 
solution, moreover, for the series expansion to make sense requires that y is small 
and hence that 0 is large, which leads to that H is large, i.e., for the present models 
the Hubble slow-roll expansion attempts to approximately describe the dynamics at 
early times for a certain solution, but which one? 
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In terms of T and X, and therefore in T and 6, the definitions in eq. 0 lead to that 
the Hubble expansion (70) can be written on the form (from now on we set k — 1) 


1 ss X‘”f(y) = cos 2 " ftf(y) where y = | N ) = ( 


(72) 


(as usual, T = T/(l — T)) where each order yields an implicit relation for a curve 
on the scalar field boundary, although it can be used to provide a series expansion 
for A" in terms of T for small T and hence for small T. 


Note that just as in the previous approximation scheme, the zeroth order Hubble 
expansion just gives X 2n = 1, which, as stated earlier, is not a good approximation 
for the attractor solution. Nevertheless, this approximation is what is sometimes 
referred to as the slow-roll Hubble approximation, see e.g. [31]. At first order, when 


f(y) = 1 + y, eq. (72) results in 


2 n 


l-X 

J2(«-1) 


COS 


2 n . 


— _ - = (-fV 

cos 2 ( n d 0 V 3 / 


(73) 


Without loss of generality, consider the neighborhood of dS + ; then the above ex¬ 
pression yields the limit (beyond the dS + point itself) 


9 



(74) 


which is the tangency condition for the center submanifold of dS + . Thus all curves 
associated with eq. (70) originate from dS±; furthermore, all curves associated with 
orders larger than zero are tangential to the attractor solutions in the small T-limit 
toward dS±. 


To improve the range and rate of convergence one can compute Pade approximants 
for the Hubble expansion for f(y), i.e., 


3H 2 » 2_(\<p) 2n [L/M]f(y) = V((j>)[L/M}f(y), 

where 

i- (i ~i)v' 


m f (y) = 
m f (y ) = 


(75) 


(76a) 

(76b) 


In Ref. [12], these approximations have been compared with the numerically com¬ 
puted attractor solutions in detail for n = 1, see Fig. 7 and Table 3 of [I2]p*| 

13 In P2J the authors used the slow-roll parameters to introduce so-called Canterbury approxi¬ 
mants in order to improve the range and accuracy of the approximations. We find that Canterbury 
approximants give more cumbersome and poorer approximations than the corresponding Pade ap¬ 
proximants when they differ (to lowest order they agree with each other). 
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Let us now turn to what is referred to as the slow-roll approximation in e.g. 

This approximation is obtained by inserting the lowest order Hubble expansion 
approximant H = s/V ((/>)/3 into 


= -2 


‘Oil 

00’ 


which yields 


-J 


(77) 


(78) 


in the present case with a monomial potential. Expressed in terms of the variables 
Ej- and X, this gives 

an 

(79) 


n 


-'—X n ~ l T, 


which when squared leads to 

E 2 


1 - X 


2 n 


COS 


2 n 


e 


X^n-l) X 2 ( n “ 1 ) COS 2 ( n_1 ) 9 


n~ 

-T 

3 


(80) 


Note that this expression coincides with the first order Hubble expansion approxi¬ 


mation given in eq. (73), and thus the slow-roll approximation is an approximation 


to the attractor solution, since it yields a curve that is tangential to the center 
manifold in the limit of small T. Thus, apart from the zeroth order approximation, 
Hubble expansion slow-roll approximations, as well as all slow-roll Hubble expansion 
based approximants, are approximations for the attractor solution. 

The slow-roll approximation can be generalized by inserting higher order Hubble 
approximants into f(y). A complication arises since different approximations can be 
obtained by multiplying expressions by various combinations of 1 ~ X 2n f(y) (they 
are all tangential to the center manifolds in the limit of small T); we here choose 
a combination that yields the above slow-roll approximation (which was chosen to 


coincide with the first order expression in (73)) when setting f(y) = 1 = g(y) below: 


n 


-X n ~ l T 


g(y ) 

\[JW 


a(y ) = f(y) - -j- 

n dy 

= i + 'i--)y+(l--)\ 2 +(i-~, r 

7i J \ nj \ n. 

n n 6 j \ n 



(81) 


Inserting the rational approximants (76) into (77) leads to the corresponding slow- 


roll approximants based on the Pade approximants constructed from the Hubble 
expansion. In Ref. [12] the corresponding curves in the state space S on the scalar 
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field subset Q m = 0 are given for n — 1 in Fig. 10 and the associated relative errors 
at the end of inflation at q = 0 are given in Table 4. Next we will compare some 
of the simplest, and thereby the most easily applicable, approximations for each of 
the above schemes. 


4.4 Approximant comparisons 


As seen from the above discussion, series expansion approximants and slow-roll 
approximation schemes lead to ambiguities, i.e., different possible types of approx¬ 
imations. This is to be contrasted with the center manifold analysis, which yield 
a systematic scheme for obtaining a unique sequence of converging approximations 
for the attractor solution. On the other hand, the strength of slow-roll approxima¬ 
tion schemes is that they can be applied to any potential and problem that satisfies 
the slow-roll conditions; it remains to be seen to what extent this can be accom¬ 
plished with other methods like center manifold analysis, an issue we will return to 
elsewhere. 


To keep the discussion reasonably short, we here restrict it to comparing only some 
of the fairly simple approximations, namely: (i) the slow-roll approximation given 
by Xq fa — |A" n_1 T, (ii) the center manifold Pade approximant 9 a — q^Tfl/l]^ 
given in eq. (66), which we will refer to as [l/lje, and 
approximant for Xq given in eq. (68), which we will refer to as [1/1] £ 
racy of the different approximations for the attractor solution, which is computed 
numerically, is shown in Figure [4| while the relative errors at q = 0 are given in 
Table HTF1 


iii) the series expansion Pade 
The accu- 


n 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

15 

20 

Slow-Roll 

14% 

3% 

19% 

33% 

44% 

54% 

62% 

69% 

75% 

79% 

92% 

97% 

me 

7% 

2% 

6% 

13% 

20% 

27% 

34% 

40% 

45% 

50% 

69% 

81% 

[T/Tjst 

14% 

19% 

22% 

24% 

25% 

26% 

27% 

27.8% 

28.3% 

29% 

30% 

31% 


Table 1: Relative Hubble errors PP at q = 0 for the slow-roll Xq a — |A" n 1 T, and 
the [1/1] 0 and [l/l]s t Pade approximants. 


All the presently discussed approximations, including the slow-roll approximation 
Xq ~ — |A n ~ x T, have quite small errors when q a — 1 , which is a reflection of 


14 


The relative errors are given in terms of the Hubble variable H by \AH\/H = 


I H num -H a 


H ,, 


pprox | 


(k 


»(1 ~ T a . 


where the subscript num stands for the numerically computed attractor 


1 \T appr o X (,l-Tnn m ) j 

solution while the subscript approx stands for the approximant. 
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(c) n = 5. (d) n = 6. 

Figure 4: Relative Hubble errors as functions of q for n = 2, 3,5, 6 for the slow- 
roll approximation (dashed curves), the center manifold Pade approximants [l/l]# 
(dash-dotted curves), and the [l/l]s t Pade approximants (space-dotted curves) 

that they are all tangential to the attractor solution when T —» 0. For larger q 
the slow-roll approximation give fairly good results for small n while it becomes 
less competitive for steep potentials with large n (recall that 4> —> ±cx) in the de 
Sitter fixed point limit), as can be expected. However, for small n the slow-roll 
approximation lead to fluctuations in accuracy when q is no longer close to — 1, 
indicating that it sometimes might be rather hard to predict how good the slow-roll 
approximation actually is going to be in a given situation. 

The above suggests that the slow-roll approximation, as well as its more accurate 
but more complicated higher expansion slow-roll based approximants, should be 
used with care when it comes to precision cosmology. It should also be pointed out 
that all approximations get worse for large n, hinting at that steep potentials might 
pose a challenge for approximative methods (although the physical viability of such 
models might be questioned in view of recent observational data). However, the 
[l/l]s t Pade approximant does not get worse so fast as the other approximations 
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and for n > 6 it is the best approximation of the three. This suggests that a good 
analytic description, pertinent in the context of precision cosmology, might require 
the development of a variety of complementary approximation methods that can be 
applied to a range of different problems; we emphasize that the present problem is 
just an illustrative stepping stone for such developments. 


5 Concluding remarks 

In this paper we have used the inherent physical scales associated with the problem 
at hand to produce new variables. In particular we obtained a regular dynamical 
system ([7]) on a compact state space S. This allowed us to give a pictorial descrip¬ 
tion of the entire solution space and moreover give proofs concerning the solutions’ 
asymptotic properties, which in turn were tied to physically important features such 
as asymptotic manifest self-similarity and breaking thereof. In addition we showed 
that the so-called attractor solution is just the unstable center manifold of a de Sit¬ 
ter state on the unphysical past boundary subset on S. We also found it convenient 
to introduce several auxiliary dynamical systems, adapted to the physical properties 
associated with asymptotic past dynamics, since this allowed us to efficiently com¬ 
pute approximations for the attractor solution by means of several complementary 
approximation schemes. The common feature of these schemes is that to lowest 
order toward the de Sitter point (apart from those approximations like the lowest 
order Hubble expansion slow-roll approximation, which just give vertical lines at 
q — — 1), from which the attractor solution originate, the approximations yield the 
tangent space of the center manifold. As a consequence they all lead to that the 
variable T along the attractor solution in the limit of small T is governed by the 
approximation 

dr 3 

This, in combination with the definition T = cH~ « and that r can be identified as 
In a, yields that when the scale factor a —>■ 0 then 

Hoc. [ln(a 0 /a)] 2 —> oo, (83) 

which is to be contrasted with the exact de Sitter solution for which H = constant. 

Finally, we would like to point out that in principle the methods we have introduced 
have a broader range of applicability than the presently studied models. For exam¬ 
ple, the use of averaging techniques should be applicable to the dynamics at late 
times for solutions that approaches a local minimum of a scalar field potential of 
the form 

V oc 0 2n ( 1 + W(4>)) where lim W — 0, (84) 

0-*- o 

where we for simplicity have translated the minimum of the potential so that it 
occurs at </> = 0, and where W obeys suitable differentiability conditions when 
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<f> —» 0. One might question the physical relevance of these potentials at very late 
times, but it should be pointed out that asymptotic dynamics can be used as pieces 
to provide approximations for global dynamics, as described in [13- 
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